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1  INTRODUCTION 


The  spherical-wave  source  scattering-matrix  description  of  electroacoustic  transducers  [1]  and 
antennas  [2]  is  ideally  suited  for  the  analysis  of  arrays  of  small  radiators  and  scatterers.  In  [1], 
it  formed  the  basis  of  the  anal3rsis  that  determined  a  simple,  expUcit,  closed-form  expression 
for  the  propagation  constants  {kd  —  I3d  diagram)  of  the  traveling  waves  on  an  infinite  linear 
periodic  array  of  isotropic  acoustic  scatterers.  And  in  [2],  it  was  applied  to  obtain  implicit 
transcendental  expressions  for  the  kd  —  fid  diagrams  of  the  traveling  waves  on  infinite  linear 
periodic  arrays  of  electric  (or  magnetic)  dipole  scatterers  oriented  perpendicular,  parallel, 
and  skew  to  the  direction  of  propagation  of  the  traveling  waves.  The  analyses  in  these  two 
previous  applications  of  the  spherical-wave  source  scattering-matrix  description  simplified 
greatly  by  assuming  the  scatterers  were  small  enough  that  only  the  fields  with  the  lowest 
order  radial  dependence  were  significant. 

In  the  present  report,  we  apply  the  spherical-wave  source  scattering-matrix  description 
to  obtain  an  implicit  transcendental  equation  for  the  propagation  constants  of  the  traveling 
waves  that  can  exist  on  an  infinite  linear  periodic  array  of  lossless  penetrable  (magnetodi¬ 
electric)  spheres.  Although  a  firamework  is  presented  for  all  orders  of  vector  spherical  waves, 
only  the  electric  and  magnetic  dipole  waves  are  included  in  the  detailed  analysis.  The  exclu¬ 
sion  of  higher-order  vector  spherical  waves  restricts  the  scope  of  this  paper  either  to  spheres 
sufficiently  small  that  only  the  dipole  scattered  fields  can  be  excited,  or  to  frequencies  for 
which  all  the  scattered  multipole  fields  are  negligible  except  the  dipole  fields.  In  all  the  par¬ 
ticular  numerical  examples  we  shall  consider,  the  magnitudes  of  the  dipole  wave  scattering 
coefficients  are  much  larger  than  those  of  the  higher-order  multipole  wave  coefficients,  and 
consequently  the  dipole  vector  spherical  waves  give  an  accurate  description  of  the  traveling 
waves  that  can  exist  on  the  array.  It  should  be  kept  in  mind,  however,  that  there  may 
well  be  traveling  waves  associated  with  the  higher-order  vector  spherical  waves  for  arrays  of 
penetrable  spheres  in  the  general  range  of  sizes  we  consider  in  some  of  our  numerical  exam¬ 
ples.  Although  three-dimensional  arrays  of  periodic  penetrable  spheres  have  been  analyzed 
in  the  past  [3],  [4],  [5],  the  present  analysis  concentrates  on  determining  kd  —  fid  diagrams 
for  the  traveling  waves  on  a  single  Unear  periodic  array  of  penetrable  spheres  with  arbitrary 
separation  and  arbitrary  values  of  real  permeabiUty  and  permittivity. 

The  report  begins  with  a  brief  review  of  the  spherical-wave  source  scattering-matrix  de¬ 
scription  of  electromagnetic  scatterers  and  then  goes  on  to  simplify  the  equations  for  spheres 
that  involve  only  electric  and  magnetic  dipole  vector  spherical  waves.  Next  these  simplified 
equations  are  applied  to  an  infinite  Unear  periodic  array  of  these  spheres.  The  resulting 
equations  are  combined,  rearranged,  and  reduced  to  a  single  transcendental  equation,  which 
is  solved  numerically,  to  obtain  kd  —  fid  diagrams  for  a  number  of  different  values  of  the 
relative  permeabiUty  and  permittivity  of  the  spheres  in  the  infinite  Unear  array.  Two-layered 
(coated)  spheres  axe  considered  as  well  as  homogeneous  spheres.  Notably,  the  kd  —  fid  di¬ 
agrams  predict  that  both  forward  and  backward  dipolar  traveling  waves  can  be  excited  on 
Unear  periodic  arrays  of  penetrable  spheres,  depending  on  the  frequency  and  the  values  of 
the  permeabiUty  and  permittivity  of  the  spheres.  In  contrast,  we  showed  in  [2]  that  only 
forward  traveling  waves  are  possible  on  Unear  periodic  arrays  of  unloaded  short-wire  dipoles. 


1 


2  SOURCE  SCATTERING-MATRIX  DESCRIPTION 
OF  A  GENERAL  SCATTERER 


Consider  a  general  scatterer  bounded  by  the  surface  S  as  pictured  in  Figure  1.  The  Maxwell 
equations  that  govern  the  time-harmonic  electromagnetic  fields  >  0)  in  the  source- 

free  free  space  between  the  surface  S  and  a  spherical  surface  Sr  of  radius  r  enclosing  the 
scatterer  and  extending  to  the  nearest  external  source  are  [6,  sec.  7.1] 

V  X  E(r)  =  ifcZoH  (la) 

V  X  H(r)  =  -ifcE(r)/Zo  (lb) 

which  imply  that 

V  •  E(r)  =  0  (2a) 

V  •  H(r)  =  0  (2b) 

where  E(r)  and  H(r)  are  the  electric  and  magnetic  fields,  Zq  is  the  free  space  impedance, 
A:  =  w/c  =  27r/A  is  the  free  space  propagation  constant,  A  is  the  free-space  wavelength,  and  c 
is  the  speed  of  light.  Eliminating  either  H(r)  or  E(r)  from  these  equations  yields  the  vector 
wave  equations 

-V  X  V  X  E(r)  -t-  fc2E(r)  =  0  (3a) 

-V  X  V  X  H(r)  +  k^H{T)  =  0  (3b) 

which  imply  ^  ^ 

V2E(r)  -I-  k^E{r)  =  0,  V  •  E(r)  =  0  (4a) 

V=*H(r)  +  Jfc=*H(r)  =  0,  V  •  H(r)  =  0.  (4b) 

Since  the  electric  field  satisfies  (4a)  in  the  source-free  free  space  region  between  the  surface 
S  and  the  spherical  surface  Sr,  the  electric  field  in  this  region  can  be  expanded  in  a  complete 
set  of  divergenceless  vector  spherical  wave  functions  [6], [7] 


EW  =  f:  E  [»£M2(r)  +  6aMS{r)  +  a<2NSW+CNSW]- 


li=l 


The  spherical- wave  expansion  of  the  magnetic  field  then  follows  from  (la) 


(5) 


OO  I 


H(r)  =  -iyoE  E 


(6) 


where  Fo  =  l/-^o  is  the  free-space  admittance.  The  vector  spherical  harmonics  MS^(r)  and 
NSS(r)  are  defined  in  Appendbc  A.  They  satisfy  the  equations 


N!S(r)  =  Iv  X  (r)  (7a) 

MS(r)  =  iv  X  N!«(r).  (Tb) 
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The  Mj^(r)  spherical  wave  functions  have  no  radial  component.  Because  of  that  the  M|^(r) 
functions  in  (5)  and  the  NS^(r)  functions  in  (6)  are  often  referred  to  as  TE  waves  while  the 
N|2(r)  functions  in  (5)  and  the  M|^(r)  functions  in  (6)  are  often  referred  to  as  TM  waves. 
The  M|j2(r)  and  N[^(r)  functions  have  the  radial  dependence  ji{kr)  where  j,  is  the  spherical 
Bessel  function  of  order  /  and  are  referred  to  as  “incoming”  waves.  The  M|^(r)  and  N|^(r) 
functions  have  the  radial  dependence  where  hj  ^  is  the  spherical  Hankel  function  of 

the  first  kind  of  order  I  and  are  referred  to  as  “outgoing”  waves.  The  expansions  in  (5)  and 
(6)  are  known  as  multipole  expansions.  The  /  =  1  terms  are  the  dipole  terms,  the  I  ^2 
terms  are  the  quadrupole  terms,  etc.  The  summations  containing  the  modal  coefiicients 
in  (5)  and  (6)  equal  the  electromagnetic  fields  produced  by  the  sources  that  reside  outside 
the  surface  Sy.  The  summations  containing  the  modal  coefficients  in  (5)  and  (6)  are  the 
outgoing  electromagnetic  fields  scattered  from  the  surface  S. 

The  “outgoing”  modal  coefficients  {b\^)  of  the  scatterer  are  related  to  the  “incoming” 
modal  coefficients  (oj^)  of  the  scatterer  by  a  linear  matrix  transformation  termed  the  “source 
scattering-matrix  equation”  [8],  [9]  for  the  spherical- wave  representation  of  the  scatterer. 
This  source  scattering-matrix  equation  can  be  written  as 


M  _ 

‘'Jm 


—  2^  2^  2^  ^lTn,l>m' 


(8) 


l'=l  «'=1 


The  coefficients  ™  source  scattering  matrix  embody  the  scattering  properties  of 

the  scatterer. 

The  spherical-wave  source  scattering  matrix  given  in  (5)— (8)  for  scatterers  differs  from 
the  classical  spherical-wave  scattering-matrix  given  in  [10,  secs.  9.18-9.24]  because  the  clas¬ 
sical  scattering  matrix  uses  radial  basis  functions  instead  of  This  small 

mathematical  change  from  to  ji  appreciably  simplifies  the  analysis  by  relating  all  the 
scattering-matrix  coefficients  directly  to  the  physical  sources. 


3  DIPOLAR  SCATTERING  FROM  LOSSLESS 
SPHERES 

We  now  specialize  our  general  source  scattering-matrix  framework  to  scattering  from  a  loss¬ 
less  penetrable  sphere  of  radius  a  with  relative  (to  free  space)  permittivity  and  relative 
permeability  fir-  The  center  of  the  sphere  is  assumed  located  at  the  origin  of  a  spherical 
polar  coordinate  system  as  shown  in  Figure  2.  The  radial  dependence  of  the  incoming  vector 
spherical  wave  functions  M2(r)  is  j,{kr)  whUe  that  of  the  incoming  wave  functions  N|J(r) 

is  ji{kr)/{kr)  for  the  radial  component  and  ^ angular  components. 
As  kr  —*0 


.  /,  \  kr-*0 

3i\kr)  ~ 


{krf 


l-3-5---(2/-Hl) 


(9a) 


4 


X 

Figure  2:  Schematic  diagram  of  a  lossless  sphere  of  radius  o  with  relative  permittivity  Cr 
and  relative  permeabihty  fir,  with  center  at  the  origin  of  a  spherical  polar  coordinate  system 
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jijkr)  kr^o  (krY  ^ 
kr  1  •  3  •  5  "  •  (2/  +  1) 


1  d 
kr  d(fcr) 


[(fcr)jj(A;r)] 


fcr->0 

r\j 


1 


/  +  1 

3  •  5  *  •  •  {21  +  1) 


{kry-'^ 


(9b) 

(9c) 


so  that  only  the  incoming  NiJj^(r)  wave  functions  are  non-zero  at  the  center  of  the  sphere. 
We  assume  that  either  the  sphere  is  sufficiently  small  that  only  the  dipole  scattered  fields  are 
excited,  or  the  frequency  is  such  that  all  the  scattered  multipole  fields  are  negligible  except 
the  dipole  fields.  Then  from  the  Mie  theory  of  scattering  from  a  sphere  [11,  sec.  9.25],  if  the 
scattered  electromagnetic  field  is  that  of  a  magnetic  dipole 


[E^c,  n»c|  ^  ^  -iFoNS] 


(10a) 


the  incident  field  that  excites  this  magnetic  dipole  is 

[E‘’“,H'n«]  =  [MSii,-iyoNS] 

whil©  if  th©  scatt©r©d  el©ctroniagn©tic  field  is  that  of  an  electric  dipole 

[E-,  H-]  =  briNfl  -iYoM^^] 
the  incident  field  that  excites  this  electric  dipole  is 

=  INS,  -iloMi'i) 


(10b) 


(10c) 


(lOd) 


m 

by 


—1,0, 1.  The  Mie  dipole  scattering  coefficients  af  and  in  (10a)  and  (10b)  are  given 


<  =  - 


6*1^=  =  — 


Hrji{^ka)[ika)ji{ka)]'  -  ji{ka)[imka)ji{mka)y 
fj,rji{'taka)[{ka)h?{ka)y  -  hY\ka)[{mka)ji{mka)]' 
tJirjiika)[{mka)ji{mka)]'  -  m^ji{mka)[ika)ji{ka)Y 
HrfYi\ka)[{mka)ji{mka)]'  —  ih^ji{mka)[{ka)h!i\ka)]' 
where  m  is  the  relative  refractive  index  m  =  Note  that  the  one-to-one  relation  of  the 

scattered  field  wave  functions  to  the  incoming  field  wave  functions,  and  hence  the  restriction 
of  the  relevant  spherical  wave  functions  to  those  with  /  =  1  and  the  same  values  of  the 
index  m  as  those  of  the  assumed  scattered  dipole  fields,  is  a  consequence  of  our  considering 
scattering  by  a  homogeneous  sphere.  Homogeneity  can  be  replaced  with  radial  symmetry 
and  the  same  conclusion  holds.  If,  however,  we  were  to  consider  a  scatterer  of  non-spherical 
shape,  or  for  that  matter  a  penetrable  sphere  with  a  non-radial  variation  of  permittivity  or 
permeability,  then  even  if  the  only  scattered  electromagnetic  field  were  that,  say,  of  (10a), 
the  exciting  incident  field  would  no  longer  be  limited  to  (10b)  and  would  in  general  include 
higher-order  wave  functions  as  well.  (Conversely,  for  these  asymmetric  scatterers,  an  incident 
dipole  field  of  the  form  (10b)  or  (lOd)  can  excite  higher-order  multipole  fields  in  addition  to 
the  dipole  fields  of  (10a)  and  (10c).) 

Useful  properties  of  the  Mie  coefficients  are  [12,  sec.  10.21] 


of 


e-2iai) 


(12a) 
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|af  1^  =  sin^tti 

(12b) 

tanai  = 

firji{mka)[{ka)ji{ka)]'  —  ji(/i:a)[(mA:a)ji(mfca)]' 
firji{raka)[{ka)yi{ka)]'  -  yi{ka)[{mka)ji{mka)]' 

(12c) 

sin  ^arg  [of]  “  |)  =  lafl 

(12d) 

arg[of|=tan 

(12e) 

f'r=-5(i-e-'“*) 

(13a) 

|6f  1^  =  sin^A 

(13b) 

tan  /?!  =  - 

firjiika)[(mka)ji{mka)]'  -  m2ji(m/!:a)[(A;a)ii(A;a)]' 
liryi{ka)[{mka)ji{inka)Y  —  ia'^ji(mka)[{ka)yi{ka)]' 

(13c) 

sm(argm-|)  =  |6n 

(13d) 

1  r  sin2)3i 

arg|M=tan  • 

(13e) 

In  (12c)  and  (13c),  y\{z)  is  the  spherical  Neumann  function.  Similar  relationships  hold  for 
all  the  Mie  coefficients  of  lossless  scatterers  [12,  sec.  10.21].^  Note  that  the  magnitude  of  the 
Mie  scattering  coefficients  is  at  most  equal  to  one.  The  values  of  ka  for  which  the  magnitude 
of  the  Mie  coefficients  a*®  and  6*“^  equals  one  correspond  to  resonances  of  the  free  magnetic 
or  electric  dipole  oscillations  of  the  sphere,  respectively  [11,  sec.  9.25].  As  examples  of 
these  resonances,  in  Figures  3,  4,  and  5  we  show  plots  of  the  magnitudes  of  the  Mie  dipole 
coefficients  a*®  and  for  a  sphere  with  Cr  =  10,  fir  =  10,  a  sphere  with  Cr  =  40,  =  !>  and 

a  sphere  with  Cr  =  10,  fir  =  1,  respectively.  (The  Mie  electric  and  magnetic  coefficients  for 
a  sphere  with  equal  relative  permittivity  and  permeability  are  identical.)  In  our  treatment 
below  of  traveling  waves  on  a  linear  periodic  array  of  spheres,  the  resonances  of  the  Mie 
coefficients  will  play  an  important  role  since  it  is  only  when  the  spheres  are  close  to  a 
resonance  that  there  is  sufficient  scattering  from  the  spheres  to  support  a  traveUng  wave. 

Let  us  now  assume  that  at  the  center  of  the  sphere  the  incoming  electric  and  magnetic 
fields  are  polarized  in  the  x  and  y  directions,  respectively.  Then  the  incoming  electromagnetic 
field  can  be  represented  as  a  linear  combination  of  the  fields  [see  (22a)  and  (22b)] 

[Er,  Hr]  =  [M^ i(r)  +  M<;>.i(r),  -iYo  (Nri(r)  +  Ni;ii(r))]  (14a) 

[Er,  Hr]  =  [n^ i(r)  -  Nj>.i(r),  -iTo  (Mri(r)  -  .  (14b) 

^For  scattering  from  homogeneous  lossless  spheres,  the  “outgoing”  spherical  mode  with  coefficient  in 
(5)-(6)  is  excited  only  by  the  corresponding  “incoming”  spherical  mode  with  coefficient  a|^  in  (5)-(6)  so  that 
all  the  in  (8)  are  zero  except  for  (s'  =  s,/'  —  =  m).  The  general  power  conservation  relations 

given  in  [2,  eq.  (33)]  then  reduce  to  Re  H-  |  =  0,  a  relationship  that  can  be  rewritten  in  the 

form  of  (12d)  and  (13d).  Also,  (5)-(6),  (8),  and  (10)  show  that  =  af'  and  =  6f ,  m  =  -1, 0, 1. 
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|Mie  coefficient]  (dB) 


8 


|Mie  coefficient!  (dB) 


The  corresponding  scattered  fields  are  then  a  linear  combination  of  the  fields 


[Bf ,  Hri  =  [MS?-i(r)  +  mS«  ,(r),  -iVo  (NSJI,(r)  +  nS«  ,(r))] 

{15a) 

|ES',Hfl=[NSJ>..{r)-NS5,,(r). 

-iy„(MS«.(r)-M«.(r))]. 

(15b) 

Because  these  fields  will  play  a  central  role  in  what  follows  it  is  convenient  to  define  a  set  of 
“combination  modes” 

NS?-W  =  NSt.(r)-Nl«,(r). 

NS?+(r)  =  NS‘,!.,(r)  +  Na,(r) 

(16a) 

Mlt(r)  =  M<;L,(r)-MS',Vi(r). 

MS%(r)  =  MS'l,(r)  +  MS«+,(r) 

(16b) 

i  =  1, 2.  These  combination  modes  are  directly  related  to  the  vector  wave  functions 
m  used  by  Stratton  [11,  p.  416]: 

n  and 

NgKr)  = 

/fnS’iW 

(17a) 

Ni^l(r)  = 

/fnS>(r) 

(17b) 

mSIw  = 

(17c) 

MitUr)  =  i. 

/f““W 

(17d) 

where  (i)  =  1  (j)  =  1  and  (i)  =  2  ^  (j)  =  3.  Explicit  expressions  for  the  outgoing 

combination  modes  are  as  follows: 


[-^  (l  +  sin(»cos#+ 1  (1  +  ^  (cos9cos«>9  -  sm^) 

(18a) 

Na(r)  =  (1  +  sin^sin#  +  1  ^  (cos0sin<^0  +  cos^^) 

MS^l(r)  =  (1  +  (shi<^^  +  cos0cos<^0) 

Mi^(r)  =  (1  +  i^os<l>9  -  cos0sin<^0). 


The  scattered  electromagnetic  fields  of  (15a)  and  (15b),  now  written  as 

(Er.Hn  =  [MS5.(r), -iVoNl^  (r)] 

and 

|Ef,Hf]  =  [Nl»(r),-iy„M»Ur)] 


{18b) 

(18c) 


(18d) 

(19a) 

(19b) 
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are  proportional  to  the  fields  radiated  by  an  infinitesimal  y  directed  magnetic  dipole  and  x 
directed  electric  dipole  respectively. 

E/Xplicit  expressions  for  the  incoming  combination  modes  are 


Nl:l(r)  -  di.«cos.^f  +  (cosecos^e-  sin^) 

(20a) 

N!;l(r)  -  ^  sin«sin^f  +  (cos«sm.^9  +  cos^)] 

(20b) 

MS|l.(r)  =  (sin^d  +  cos  0  cos  0^) 

(20c) 

M$^l(r)  =  7:\  -ji{kr)  (cos<f>e  -  cos0sin(?i^) 

’  2  V  ^  ' 

(20d) 

cos^ 

J.W  -  ^2  ^ 

(21a) 

ji(z)  sinz  cos  2; 

z 

(21b) 

1  d  r  .  ,  M  COS2:  sin^  /  1  \ 

jgjfe.Wl-  ^ 

(21c) 

Limiting  expressions  for  the  incoming  combination  modes  bs  kr  0  are 

(22a) 

(22b) 

and  ,,,  .  „ 

Mi;i(r)'"=°0(fcr). 

(22c) 

Prom  the  Mie  theory  of  scattering  from  a  sphere  (lOa)-(lOd),  the  source  scattering-matrix 
equations  for  the  combination  modes  can  be  written  as 

e=s!;i£u:i 

(23a) 

M«=siSU« 

(23b) 

where  Cb^\.  is  the  coefiicient  of  the  incoming  electromagnetic  field  ii^Ni^^.(r)],  ^1,4. 

is  the  coefficient  of  the  corresponding  outgoing  electromagnetic  field  [MS^{.(r),  -iFoMi,+(r)], 
is  the  coefiicient  of  the  incoming  electromagnetic  field  [NSJI(r),-iyoMS5l(r)],^6i%  is 
the  coefiicient  of  the  corresponding  outgoing  electromagnetic  field  [Nl^L(r),  -iloMi,i(r)], 
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and  the  scattering  coefficients  and  are  equal  to  the  Mie  coefficients  given  by 

(11a)  and  (11b) 

=  aT  (24a) 

=  (24b) 

Note  firom  Footnote  1  that  5i^+fi,+  =  and  rn  =  ±l. 

It  will  be  convenient  to  normalize  the  scattering-matrix  equations.  As  r  — 0  the  incoming 
electric  field  E®(r)  =  ui^lNi^l(r)  reduces  to 

E“(r)’-3^o®-i=*  (25a) 

y/oTT 

and  the  incoming  magnetic  field  H°(r)  =  ai|+(— il^)Ni|^(r)  reduces  to 

H“(r)  a?i(-iFo);^y.  (25b) 


The  incident  electric  and  magnetic  dipole  spherical  waves  are  the  only  multipoles  with 
nonzero  fields  at  the  center  of  the  sphere  (r  =  0).  Let 


Then 


a?_  = 

1 


®1,+  -  ~i 


Now  let  6_  be  the  coefficient  of  - —  in  the  transverse  component  of  the  outgoing  electric 

kr 

glfcr 

field  6i^lNi^^I(r)  [see  (18a)]  and  6+  be  the  coefficient  of  in  the  outgoing  electric  field 


(28a) 

(28b) 

Then  the  scattering  matrix  equations  (23a)  and  (23b)  become 

b-  =  -i|  S'Plfl  E° 

(29a) 

b+  -  I2  *l.+;l,+  • 

(29b) 

Finally,  normalize  the  scattering  coefficients  by  letting 

S-  =  S.%?-  =  hi 

(30a) 

13 


(30b) 


C  _  cO).(l) _ i- 

0-1-  —  —1—  Oj  _j_.i 

Then  the  normalized  form  of  the  scattering  matrix  equations  is  simply 

b-=-S.E^ 

b+  =  5+ 

To 

Since  these  equations  form  the  basis  of  our  treatment  of  dipolar  traveling  waves  on  an 
infinite  linear  periodic  array  of  penetrable  spheres  it  is  worthwhile  to  elucidate  their  meaning. 
Equation  (31a)  says  that  if  the  electric  field  incident  on  a  penetrable  sphere  has  the  value 
E®  X  at  r  =  0,  then  the  corresponding  scattered  electric  and  magnetic  dipole  fields  are  [see 

(18a)  and  (18c)] 


(31a) 

(31b) 


=  b- 


nXkr 


kr 


(' + s)  + (' + F  -  (]^) 

/  i  \ 

cos6cos<f)<f>). 


(32a) 

(32b) 


Equation  (31b)  says  that  if  the  magnetic  field  incident  on  a  penetrable  sphere  has  the  value 
y  at  r  =  0,  then  the  corresponding  scattered  electric  and  magnetic  dipole  fields  are  [see 
(18d)  and  (18b)] 


=  64 


( 


=  -iYob+ 


rSkr 


kr 


(cos  <f>0  —  cos  $  sin  <f)^)  (33a) 

krj 

A  ^i  +  ±^sin0sin(^f +  i  ^1  +  ^  -  (cos0sin«^0  +  cos0^) 


(3 


lb) 


The  scattered  fields  in  (32)  and  (33)  are  those  of  an  a:  directed  electric  dipole  and  y  directed 
magnetic  dipole,  respectively.  In  treating  on-axis  scattering  from  an  infinite  linear  periodic 
array  of  penetrable  spheres  in  the  next  section  it  will  be  useful  to  have  expressions  for  the 
special  cases  of  equations  (32)  and  (33)  when  9  =  0  and  9  =  ir,  that  is,  on  the  axis  of 
the  array.  Prom  the  standard  expressions  relating  the  unit  vectors  in  spherical  polar  and 
cartesian  coordinates  we  obtain  from  (32) 


(34a) 

(34b) 

and  from  (33)  we  obtain 

(35a) 

(35b) 
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4  ON-AXIS  DIPOLAR  SCATTERING  EQUATIONS 
FOR  AN  INFINITE  LINEAR  PERIODIC  ARRAY 
OF  LOSSLESS  PENETRABLE  SPHERES 


The  on-axis  dipole  scattering-matrix  equations  for  a  infinite  linear  periodic  array  of  lossless 
penetrable  spheres  with  radius  a,  relative  permittivity  €r  and  relative  permeability  fir  can 
now  be  easily  obtained.  These  equations  are  the  basis  for  the  treatment  of  dipolar  traveling 
waves  on  an  infinite  linear  periodic  array  of  lossless  penetrable  spheres  given  in  Section  5. 
We  denote  the  separation  between  the  centers  of  adjacent  spheres  by  d,  take  the  z  axis  to  be 
the  axis  of  the  array,  and  assume  an  excitation  of  the  array  with  an  x  directed  electric  field 
and  a  y  directed  magnetic  field  (see  Figure  6).  The  x  directed  electric  field  at  the  center 
of  the  nth  sphere  resulting  firom  the  electric  fields  scattered  by  all  the  other  spheres  in  the 
array  is  then,  from  (34a)  and  (35a), 

00  giMi-nt  /  i  1  \ 

~  kd\j  —  n\  \  ^  kd\j  —  n\  (fcd)2|j  — n^y 


i=-oo 


n-1  ^ikd{n-j)  / 

Since  6_,n  =  S-E^ 


H- 


•  \  00  giA:d(j-n)  /  j  \ 

kd{n  -  j) )  kd{j  -  n)  ^  kd{j  -  n)  j  ' 


r  00  gifedlj-«l  /  i  1 

i  '’-■‘kd\3-n\  kd\j-n\  ~  (kd)^\j-np) 

j^n 

n-l  pifcd(n-i)  /  i  \  /  i  \l 

^'^'^kd{n-j)  kd{n-j))  “  ^-^’^kd{j-n)  kd{j  -  n)  j  J  ‘ 

Similarly  the  y  directed  magnetic  field  at  the  center  of  the  nth  sphere  resulting  from  the 
magnetic  fields  scattered  by  all  the  other  spheres  in  the  array  is,  fi-om  (34b)  and  (35b), 


ffOn  00  pife<i|j-nl  / 


n-l 

+ 

jsr—OO 


Qikd(n-j) 

kd{n  -  j) 


j=-oo 

j^n 


1  + 


kd\j  —  n|  {kdy\j  —  n| 


i  \  00  giferfO-n)  /  i  \ 

^  kd{n  -  j) )  kd{j  -  n)  V  ^  kd{j  -  n)  j  ‘ 


ROn 

Since  6+,„  =  S+^ 

iO 


<>+,«  =  s+  E  <> 


^  j=—oo 


(■* 


Qikd\j-n\ 

kd\j  —  n]  '  kd\j  —  n|  {kdy\j 


3  -  «l^) 


+ 


n-l  Jkd(n-j)  (  :  \  oo  Akd(j-n)  (  i  \  1 

Kn-J)  kd{n-j))  (‘  +  kdU-n)]  I 


15 
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Figure  6:  Schematic  diagram  of  an  infinite  linear  periodic  array  of  lossless,  penetrable  spheres 
of  radius  a  with  relative  permittivity  €r,  relative  permeability  fir,  and  separation  distance  d. 
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5  DIPOLAR  TRAVELING  WAVES  ON  AN  INFINITE 
LINEAR  PERIODIC  ARRAY  OF  LOSSLESS  PEN¬ 
ETRABLE  SPHERES 


We  now  come  to  the  central  objective  of  this  report  —  the  derivation  of  an  expression  for 
the  propagation  constant  of  a  dipolar  traveling  wave  supported  by  an  infinite  linear  periodic 
array  of  lossless  penetrable  spheres.  If  the  array  supports  a  dipolar  traveling  wave  with  a 
real  propagation  constant  jS  to  be  determined,  then  the  b-j  in  (37)  and  (39)  are  identical 
except  for  a  phase  shift  given  by 

(40a) 
(40b) 


b-j  =  0<  fid  <  IT 

and  similarly 

b+j  =  6+,oe‘^^^  0  <  fid  <n. 
This  allows  (37)  and  (39)  to  be  written  as 

jVn 


kd\j  -  n\  {kd)'^\j 


j-np) 


+V±[  ”f''  fn.  i  \ 


oo  f3}kd(j-n)  /  ;  \  -i 


and 


1  =  5- 


1,^0  kd\j-n\{ 


1  + 


kd\j  -  nl  {kdy\j 


j-np) 


jVn 

^  bo,+ 1- ji'oo  kd{n-j)\  kd{n  -j)) 

OO  Akd(J-n)  /  \  \  1 


Since  the  values  of  the  summations  are  independent  of  the  value  of  n,  we  can  set  n  =  0. 
Multiplying  by  {kdtf  we  then  obtain  the  equations 

00  /  ^  Q-i{p-k)d3^ 


M»  =  S-[e( 

^0,+ 


+ 


K- 


oo  /pi(fc+/3)4j  _  (.-K0-k)^ 

ME 


i=i 


(43) 


and 


(kd)- 


00  /f,\(k+fi)dj  ^  Q-i(0-k)cb 


)( 


,,  ,.o  ikd  1 
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^0, 


Hkd)t{ 


oo  /gi(fc4-i9)c&‘  _  g-i(/3“fc)4/ 

,“i  \  3 


Now  let 


P 


0,+ 


^0.- 


(44) 


(45) 


(Recall  —  see  (23a,  23b),  (15a,  15b),  and  (40a,  40b)  —  that  the  electromagnetic  field  scat¬ 
tered  from  each  of  the  spheres  in  the  array  is,  apart  from  the  phase  factor  of  the  traveling 
wave,  a  linear  combination  of  an  electromagnetic  field  proportional  to  the  field  of  an  infini¬ 
tesimal  X  directed  electric  dipole  with  a  normalized  coefficient  6_,o,  and  an  electromagnetic 
field  proportional  to  the  field  of  an  infinitesimal  y  directed  magnetic  dipole  with  normalized 
coefficient  6+,o,  so  that  p  is  the  ratio  of  these  two  normalized  scattered  field  coefficients.) 
Also  let 


Si 


00  /^iik+M  _|. 


=e( 

j=l  \ 


and 


=  A:df:(- 

i=i  \ 


.i(fc+/3)4»  _  Q-K0-k)<lli 


Using  the  summation  formulas  [13,  sec.  1.441,  eqs.  1  and  2;  sec.  1.443,  eqs.  3  and  5] 

1 


<11^-— s  COSJffl  1 1 

Pi  j  ”2  2(l-cosa) 


=  —  In  [2  sin(o/2)] ,  0  <  o  <  27r 


,  sin  ja  TT  —  a 


^  i 

3=1  J 


,  0  <  a  <  27r 


3=1  3 

and  the  approximations  (see  Appendix  B) 
sin  jo 


~co^^^_7ro  ^  0<o<27r 

4^  6  2  4 

^  sin  70  TT^o  -nc?  ^ 

— 4-  - h  Tpr,  0  <  o  <  27r 

-3  fi  4  ^  12’ 


(46a) 

(46b) 

(47a) 

(47b) 

(47c) 

(47d) 


j,(a)fa_o.i381sino-|-0.03212sin2o-0.9653oln(o/7r),  0  <  o  <  tt  (48a) 

3=1  3 


F{a)  =  -F{2'n  -a),  tt  <  o  <  27r 


(48b) 


and 

OO 


=  G^a)  RJ  1.3328-0.1424  cos  a-|-0.01094  cos  2a+0A902a^  ln(o/7r)-0.2417o2,  q  <  o  <  tt 

(48c) 

G(o)  =  G(27r  -  o),  7r<o<27r  (48d) 
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the  summations  Ei  and  E2  can  be  expressed  in  closed  form  as 


4 


i 


Si  =  —{kdf  In  [2  (cos  kd  —  cos  /?d)]  —  (fed)  [F{kd  +  fid)  —  F{Pd  —  kd)] 


-  [G{kd  +  pd)  +  G{pd  -  kd)]  -  i-(fcd)' 


(49) 


and 

S2  =  -(fed)* 


Insin^^^^^  -Insin^^^^^j  -(fcd)[F(A:d  +  |0d)  +  F(/3d-fcd)]. 


(50) 

(51a) 


The  on-axis  scattering  equations  (43)  and  (44)  can  then  be  written  as 

(kdf  =  S-  (El  -I-PE2) 

and  .  . 

(fcd)^  =  5+  f  El  -h  ^£2  j  . 

Equations  (51a,  51b)  are  a  pair  of  simultaneous  equations  determining  the  two  unknowns 
p  and  Pd.  We  can  solve  for  p  in  each  of  these  two  equations  and  equate  the  resulting 
expressions,  thus  obtaining 


(51b) 


(kd)^  -  5_Ei 
5_E2 


5+S2 


(fcd)3-5+Ei' 


(52) 


It  is  important  to  note  that  the  quantity  p  given  by  either  side  of  (52)  is  real.  For,  using  the 
expression  for  p  given  by  the  LHS  of  (52),  the  fact  that  the  imaginary  part  of  Ei  is  --(fed)® 
[see  (49)]  and  that  E2  is  real  [see  (50)],  the  definition  of  S-  in  (30a),  and  the  property  (13d) 
of  the  Mie  coefiicient  it  follows  that 


E2  Im[p]  =  (fcd)®  Im 


S- 


-  Im[Ei]  =  {kdf 


(53) 


The  same  conclusion  can  be  obtained  similarly  firom  considering  the  RHS  of  (52).  Equation 
(52)  can  then  be  easily  solved  nmnerically  for  pd  given  kd. 

It  is  worth  noting  as  a  check  on  the  consistency  of  what  we  have  done  that  if  we  divide 
(51a)  by  S-  and  then  let  S-  =  |5_|  exp(iV»-)  and  equate  real  and  imaginary  parts  with  p 
real,  we  obtain 

^*‘^^%in^_  =  -|(jfcd)3  (54) 


|5-| 


or 


(65) 


a  result  that  is  consistent  with  (30a)  and  (13d).  A  similar  check  can  be  performed  with  (51b) 
to  produce  liS+l  =  (3/2)  sin  We  note  that  (55)  is  identical  with  the  equation  obtained 
for  the  magnitude  of  the  effective  scattering  coefficient  of  a  small,  reciprocal,  lossless  electric 
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dipole  in  [2].  As  a  further  check  on  the  validity  of  (52),  we  note  that  when  S+,  the  normalized 
magnetic  dipole  scattering  coefficient,  is  equal  to  zero,  then  (52)  reduces  to 

(kdf  -  5_Ei  =  0.  (56) 

This  equation  is  identical  to  equation  (72)  of  [2]  that  determines  the  k  ^  diagram  for  the 
traveling  waves  on  a  linear  array  of  electric  dipoles  perpendicular  to  the  array  axis. 

It  is  easy  to  show  that  if  fid  is  a  solution  of  (52)  then  so  is  27r  —  ^d  (Si  is  unchanged 
and  E2  is  reversed  in  sign  if  2'K  —  Pd  is  substituted  for  pd)  .  The  significance  of  this  is  that 
if  one  had  a  fast  wave,  for  example,  tt  <  /3d  <  27r  then,  as  shown  in  general  by  Yaghjian 
[1],  P'  with  0  <  /J'd  =  27r  -  ^d  <  TT  is  the  propagation  constant  of  a  slow  traveling  wave  in 
the  negative  ^  direction,  exp(-i)9'2)  .  Therefore,  it  is  unnecessary  to  consider  values  of  pd 
greater  than  tt  in  the  kd  -  pd  diagrams  of  traveling  waves. 


6  COATED  SPHERES 


The  analysis  that  has  been  performed  in  the  previous  three  sections  has  assumed  that  the  rel¬ 
ative  permittivity  and  permeability  of  the  spheres  is  constant  throughout  the  entire  sphere. 
The  analysis,  however,  carries  over  almost  identically  if  the  relative  permittivity  and  perme¬ 
ability  of  the  spheres  is  allowed  to  vary  from  one  spherical  shell  to  another.  As  one  of  the 
simplest  examples,  assume  that  the  relative  permittivity  and  permeability  of  each  sphere  is 
equal  to  Cri  and  fjLri,  respectively,  for  0  <  p  <  a,  and  is  equal  to  er2  and  Pr2»  respectively, 
for  a  <  p  <  5  where  b  is  the  radius  of  the  sphere.  Let  mi  and  m2  be  the  relative  refractive 
indices  mi  =  and  m2  =  respectively,  for  the  two  regions  of  the  sphere. 

Then  [14,  sect.  8.1]  the  scattering  coefficients  of  the  coated  sphere,  and  b^,  analogous  to 
the  Mie  coefficients  af  and  6f  for  the  homogeneous  sphere  given  by  (11a)  and  (11b)  are  as 
follows: 


Jl(P2)  {[ni2P2jl(ni2P2)]"  +  A  [m2p2!/l(m2P2)]'}  ~  P2r  [P2jl(P2)]'  \jl(.^2P2)  +  Ayn(gl2P2)] 

h^l\p2)  {[ni2p2jl(in2P2)]"  +  A  [m2P2yi(ni2P2)]'}  -  IJ-2t  [p2hf^(P2)]  bl(“2P2)  +  Ay„(m2P2)] 

(57) 

Pr2jl(P2)  {[m2P2ji(m2P2)]^  +  B  [m2P2yi(m2P2)]^}  -  nxj  \P23\{P2)\  b'l(”^2P2)  +  %n(m2p2)] 

Pr2hP(p2)  {[m2P2ji(ni2P2)]'  +  B  [m2P2yi(m2p2)]'}  -  mi  [p2h^i\p2)]  l7i(m2P2)  +  5y„(m2p2)] 

(58) 


where  pi  =  fca,  p2  =  kb,  and 

_  Priji(miPi)[(m2Pi)ji(m2Pi]^  —  Pr2ji(m2Pi)[niipiji(mipi)]^ 

~  Priji(miPi)[(m2Pi)yi(m2Pi)]'  -  Pr2yi(m2Pi)[(mipi)ji(miPi)]' 

_  Primjii  (m2Pi )  [(mipi  )ji  (mipi)]^  -  Pr2miii  (mipi)  [(m2Pi  )ji  (m2Pi)]^ 
”  Primiyi(m2Pi)[(niipi)yi(mipi)]'  -  Pr2mfii(mipi)[(m2Pi)yi(m2Pi)]' 

The  equations  corresponding  to  (12a-12e)  and  (13a-13e)  are  as  follows: 


(59) 

(60) 


®lc 


(61a) 
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and 


where 


and 


|af^|^  =  sin^aic 

(61b) 

^  fi2rAa*  —  aA' 

IJi2rAcf-cA' 

(61c) 

sin  ^arg  - 

i)  =  |“SI 

(61d) 

arg  [af^]  =  tan”^ 

r  sin2aic 

(61e) 

L— 2sin^aic. 

-  e-2‘^‘-=) 

(62a) 

|6i“/  =  smVu 

(62b) 

^  „  IMtrB'a  -  n^Ba' 

/i2rjB'c  —  mlBcf 

(62c) 

sin  ^arg  [6J3  - 

|)=ra 

(62d) 

arg  [*ic]  =  tan"^ 

r  sin2/?ic 

(62e) 

[-2sin^/3ic 

A  =  ji{m2p2)  +  Ayi{m2P2) 

(63a) 

M  =  [(m2P2)ji(m2/!)2)]'  +  ^[(m2P2)yi(m2p2)]'  (63b) 

a=ii(p2) 

(63c) 

0'  =  \p2ji{p2Y 

(63d) 

c  =  yi{p2) 

(63e) 

c'  =  \p2y\{p2) 

(63f) 

B  =  ji{Ta2p2)  +  Byi(m2p2) 

(63g) 

B'  =  [(m2/92)ji(m2P2)]'  +  fi[(m2/02)yi(m2/!>2)]' .  (63h) 

The  central  result  of  Section  5,  equation  (52)  that  implicitly  gives  the  propagation  con¬ 
stant  of  a  dipolar  traveling  wave  on  an  array  of  homogeneous  spheres  in  terms  of  the  free-space 
wave  number  and  the  normalized  scattering  coefficients  of  the  spheres,  then  holds  for  coated 
spheres  as  well  if  and  0^^  are  substituted  for  and  a",  respectively,  in  (30a)  and  (30b). 
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7  NUMERICAL  RESULTS 

Computer  calculations  have  been  performed  to  illustrate  the  theory  of  dipolar  traveling  waves 
on  linear  periodic  arrays  of  lossless  spheres  that  we  have  presented.  The  computations  fall 
into  two  main  categories.  In  the  first  set  of  calculations  the  electrical  radius  of  the  spheres,  ka, 
is  held  constant  and  the  electrical  inter-sphere  separation  distance,  kd,  is  varied  to  obtain  13d 
where  /?  is  the  propagation  constant  of  the  traveling  wave.  In  the  second  set  of  computations 
the  ratio  of  the  sphere  radius  to  the  inter-sphere  separation  distance,  a/d,  is  held  constant, 
and  kd  is  varied  to  obtain  /3d.  In  both  sets  of  calculations  we  present  the  kd  -  13d  diagram 
showing  a  plot  of  kd  versus  fid.  The  second  set  of  calculations  corresponds  to  a  laboratory 
measurement  procedure  in  which  frequency  is  varied  and  so  is  much  more  practicable  than 
the  first  set  which  corresponds  to  a  fixed  frequency  in  a  laboratory  naeasurement  procedure. 

To  begin  with  we  show  some  kd — fid  diagrams  for  a  linear  periodic  array  of  spheres  with 
_  iQ  for  ka  equal  to  several  different  resonances  of  the  Mie  scattering  coefficients. 
Figure  7  shows  the  kd  -  fid  diagrams  for  ka  =  0.405, 0.693,  and  0.998.  The  smallest  value  of 
kd  in  these  diagrams  must,  of  course,  be  greater  than  2ka,  otherwise  the  spheres  wiU  overlap. 
We  observe  that  the  curves  end  when  fid  =  kd,  that  is,  when  the  traveling  wave  propagation 
constant  fi  equals  the  fi:ee-space  wavenumber  k.  The  importance  of  having  fco  equal  to,  or 
close  to,  a  resonance  of  the  Mie  coefficients  and  (a*®  =  when  Cr  =  fir)  is  intuitively 
clear  since  a  traveling  wave  with  fi  significantly  greater  than  k  cannot  be  excited  without  a 
high  degree  of  scattering  coupfing  between  the  spheres  composing  the  array.  For  small  values 
of  the  scattering  coefficients  such  coupling  cannot  exist.  IVaveling  waves  with  /3  >  k  on  a 
linear  periodic  array  of  lossless  spheres  thus  exist  only  within  a  fairly  narrow  ka  “window  . 

Next  in  Figure  8  we  show  the  kd  —  fid  diagram  for  a  linear  periodic  array  of  spheres 
with  €r  =  40,  Hr  =  1  and  ka  =  0.480,  very  close  to  the  resonance  of  the  of  Mie  coefficient 
at  feo  =  0.487.  The  shape  of  the  kd  -  fid  diagram  is  somewhat  the  reverse  of  the  kd  -  fid 
diagram  of  Figure  7  for  e,  =  10,  Hr  =  10,  and  ka  =  0.405.  In  Figure  9  we  show  a  detail  of 
Figure  8  for  the  region  of  1.060  <  kd  <  1.100.  What  is  particularly  interesting  here  is  the 
region  for  1.060  <  kd  <  1.082  where  there  are  two  values  of  the  traveling  wave  propagation 
constant  for  each  inter-sphere  separation  distance. 

Figure  10  demonstrates  the  remarkable  sensitivity  of  the  kd  -  fid  diagram  of  Figure  8 
to  the  size  of  the  spheres.  As  ka  increases  only  the  small  amount  from  0.480  to  0.4844 
the  kd  -  fid  diagram  shrinks  to  a  small  curve  in  the  vicinity  of  fid  =  3.0,  kd  =  1.8.  For 
ka  =  0.4845  no  traveling  wave  can  be  supported  by  the  linear  array  of  spheres. 

Another  example  of  the  possibility  of  there  being  two  different  traveling  waves  that  can 
be  supported  by  the  same  linear  array  of  lossless  spheres  is  shown  in  the  kd  —  fid  diagram  of 
Figure  11  for  c,.  =  10,  Hr  =  1,  and  ka  =  1.1.  (The  value  of  ka  =  1.1  is  approximately  halfway 
between  the  resonance  of  the  Mie  scattering  coefficient  af  at  0.951  and  the  resonance  of  6f 
at  1.251 .)  For  2.337  <kd<  2.346  there  are  two  values  of  fid  corresponding  to  each  value  of 

the  electrical  inter-sphere  separation  distance  kd. 

We  now  turn  to  the  second  group  of  calculations,  those  for  which  the  ratio  of  the  sphere 
radius  to  the  inter-sphere  separation  distance,  a/d,  is  held  constant,  and  kd  is  varied  to  obtain 
fid.  As  noted  above,  these  calculations  correspond  to  the  practical  laboratory  measurement 
procedure  of  setting  up  an  array  of  spheres  and  then  investigating  the  behavior  of  traveling 
waves  on  the  array  as  the  frequency  is  varied.  Figure  12  shows  the  kd— fid  diagram  obtained 
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Figure  7:  kd  —  Hifl.gr a.ms  for  traveling  waves  on  an  infinite  linear  periodic  array  of  lossless, 

penetrable  spheres,  =  10,  fir  =  10. 
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Figure  10:  Sensitivity  of  kd  —  fid  diagrams  for  traveling  waves  on  an  infinite  linear  periodic 
array  of  spheres  with  =  40,  Hr  —  1>  varying  from  0.4800  to  0.4844. 


Figure  11:  kd  —  fid  diagram  for  a  traveling  wave  on  an  infinite  linear  periodic  array  of 
spheres  with  Cr  =  10,  fir  =  1,  and  ka  =  1.1. 
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Figure  12:  kd 
spheres  with  er 
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diagram  for  a  traveling  wave  on  an  infinite  linear  periodic  array  of 
10,  fir  =  10,  and  a/d  =  0.45. 
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for  a  linear  array  of  lossless  homogeneous  spheres  with  er  =  10, /ir  =  10,  and  a/d  =  0.45. 
(The  dotted  line  shows  /3d  =  kd.)  We  note  that  there  are  four  sections  of  the  fed  -  Pd 
diagram,  each  section  corresponding  to  a  narrow  window  of  ka  in  the  vicinity  of  one  of  the 
resonances  of  the  Mie  scattering  coefficients  of  and  at  ka  =  0.405, 0.693, 0.988,  and 
1.299;  see  Figure  3.  Each  section  begins  and  ends  when  fid  fn  kd.  It  is  of  considerable 
interest  that  each  of  the  four  sections  consists  of  two  branches,  the  lower  branch  with  /3d 
increasing  from  fed  to  as  fed  increases,  and  the  upper  branch  with  ydd  decreasing  from  tt 
to  fed  as  fed  continues  to  increase.  The  phase  velocity  of  the  traveling  wave  is  positive  on 
both  branches  while  the  group  velocity  (dk/dfi)  of  the  traveUng  wave  is  positive  on  the  lower 
branch  and  negative  on  the  upper  branch.  Hence  the  traveling  waves  corresponding  to  the 
upper  branches  of  the  four  sections  of  fed  —  Pd  diagram  are  “backward”  traveling  waves. 

Our  understanding  of  the  fed  -  Pd  diagram  of  Figure  12  can  be  increased  by  extending 
the  dia.gra.Tin  to  include  values  of  Pd  between  tt  and  27r.  In  Figure  13  we  show  the  lowest 
of  the  four  sections  of  the  fed  —  Pd  of  Figure  12  extended  to  pd  >  tt.  (The  dotted  line  is, 
as  in  Figure  12,  the  plot  of  /3d  =  fed.)  As  we  have  commented  at  the  end  of  Section  5,  if 
TT  <  )3d  <  27r  then  p  is  the  propagation  constant  of  a  fast  traveling  wave  in  the  positive 
2  direction,  exp(i/32:),  and  P'  with  Q  <  P'd  -  2tt  -  Pd  <  it  Ss  the  propagation  constant  of 
the  corresponding  slow  traveling  wave  in  the  negative  z  direction,  exp(— i/3'z)  .  The  dashed 
portions  of  the  plots  corresponding  to  /0d  >  tt  are  thus  simply  mirror  images  in  the  line 
/3d  =  TT  of  the  solid  portions  of  the  plots  for  0  <  pd  <  ir.  In  general  there  is  thus  no 
need  to  show  the  fed  —  pd  diagram  for  values  of  Pd  >  n  Here,  however,  what  we  can  see 
from  Figure  13  is  that  by  extending  the  fed  —  /3d  diagram  to  include  values  of  /3d  >  tt  we 
obtain  two  curves,  one  consisting  of  the  lower  solid  segment  and  the  upper  dashed  segment, 
and  the  other  consisting  of  the  upper  solid  segment  and  the  lower  dashed  segment.  Either 
curve  can  be  used  to  obtain  the  other  by  making  use  of  the  mirror  symmetry.  The  slopes 
of  these  two  curves  are  continuous  at  Pd  =  ir.  The  upper  solid  segment  corresponds  to  a 
backward  traveling  wave  since  the  group  velocity,  proportional  to  the  slope  of  the  curve,  has 
the  opposite  sign  of  the  direction  of  propagation  of  the  slow  wave.  The  group  velocity  of 
the  slow  traveling  wave  in  the  positive  z  direction  whose  fed  —  Pd  diagram  is  the  lower  solid 
segment  is  very  small  in  the  interval  between  ka  =  0.3925  and  0.3928. 

As  another  example  of  the  fed  —  Pd  diagram  obtained  when  the  ratio  of  the  sphere  radius 
to  the  inter-sphere  separation  distance  is  held  constant  and  the  frequency  varied.  Figure  14 
shows  the  fed  —  Pd  diagram  corresponding  to  Figure  12  when  =  20  =  /Xr  =  20  instead  of 
10;  o/d  =  0.45  as  before.  There  are  now  nine  sections  of  the  kd—pd  diagram  instead  of  four, 
corresponding  to  the  more  densely  spaced  resonances  of  the  Mie  scattering  coefficients  of 
af  and  6f  at  ka  =  0.214,0.366,0.516,0.665,0.816,0.968,1.122,1.279,  and  1.437.  Although 
otherwise  the  fed  —  /3d  diagrams  appear  to  be  similar,  there  is  an  interesting  new  feature  of 
the  two  lower  branches  of  the  kd—pd  diagram  of  Figure  14  not  present  in  any  of  the  branches 
of  the  kd  —  pd  diagram  of  Figure  12  or  the  higher  branches  of  Figure  14.  This  feature  can  be 
seen  clearly  in  Figure  15  which  shows  in  detail  the  behavior  of  the  lowest  branch  of  Figure  14. 
We  note  that  as  the  branch  is  traced  out  from  its  beginning  at  ka  =  0.2098  and  kd  =  0.4662 
to  its  end  at  ka  =  0.2166  and  fed  =  0.4813  that  ka  and  fed  are  not  monotonic  increasing  as 
they  are  for  all  the  branches  of  the  fed  —  pd  diagram  of  Figure  12  and  all  but  the  two  lowest 
branches  of  Figure  14.  Instead,  ka  ffist  increases  from  0.2098  to  0.2140,  then  decreases  to 
0.2107,  and  then  increases  monotonically  to  0.2166.  As  a  consequence  of  this  behavior,  in 
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Figure  13:  kd  —  fid  HingTAnn  for  a  traveling  wave  on  an  infinite  linear  periodic  array  of 
spheres  with  =  10,  /itr  =  10,  and  a/d  =  0.45. 
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Figure  14:  kd  —  fid  diagram  for  a  traveling  wave  on  an  infinite  linear  periodic  array  of 
spheres  with  =  20,/ir  =  20,  and  a/d  =  0.45. 
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Figure  15;  Lowest  branch  of  the  kd  —  fid  diagram  for  a  traveling  wave  on  an  infinite  linear 
periodic  array  of  spheres  with  Cr  =  20,  fXr  —  20,  and  a/d  =  0.45. 


the  interval  of  the  electric  spacing  of  the  spheres,  kd,  from  0.4682  to  0.4756,  there  is  not  just 
one  corresponding  value  of  the  traveling  wave  propagation  constant  Pd  but  three.  In  other 
words,  the  array  can  support  three  distinct  traveling  waves  for  kd  in  this  interval.  Also,  note 
that  this  branch  of  the  kd-pd  diagram  displays  negative  group  velocity  for  ka  between  0.2140 
and  0.2107,  and  for  ka  between  0.2110  and  0.2166,  and  positive  group  velocity  elsewhere. 
Similar  behavior,  though  less  pronounced,  is  exhibited  by  the  next-to-lowest  branch  of  the 
kd  —  Pd  diagram  of  Figure  14,  as  shown  in  Figure  16. 

The  kd  —  Pd  curves  shown  in  Figiures  12-16  have  been  for  arrays  of  lossless  magnetodi¬ 
electric  spheres.  It  is  also  of  interest  to  show  some  examples  of  kd  —  pd  curves  for  linear 
arrays  of  lossless  dielectric  spheres  =  !)•  hi  Figure  17  we  show  the  kd  —  Pd  curve  for  an 
array  of  dielectric  spheres  with  =  10  and  ajd  —  0.4.  For  small  kd  and  ka  there  is  very 
little  scattering  between  the  spheres  and  the  kd  —  Pd  curve  very  closely  follows  the  kd  —  Pd 
line,  mfianing  that  the  traveling  wave  is  very  loosely  bound  to  the  array.  When  ka  starts  to 
approach  the  first  resonance  of  the  first  magnetic  Mie  coefficient  at  ka  —  0.951  (see  Figure 
5),  scattering  between  the  spheres  increases  and  a  slow  traveling  wave  can  be  supported  by 
the  array.  The  upper  branch  of  the  kd  —  pd  curve  corresponds  to  values  of  ka  for  which 
both  the  first  magnetic  and  first  electric  Mie  coefficients  are  about  4  dB  below  their  first 
resonance  peaks.  In  Figure  18  the  kd  —  pd  curve  is  shown  for  a  linear  array  of  the  same 
dielectric  spheres,  =  10,  but  with  a/d  =  0.3.  The  lower  branch  of  the  kd  —  Pd  curve  is 
similar  in  shape  to  that  of  the  lower  branch  of  the  kd  —  Pd  curve  in  Figure  17.  The  veiry 
small  upper  branch  has  a  negative  slope  and  corresponds  to  a  backward  traveling  wave. 


8  CONCLUDING  REMARKS 

In  this  report  we  have  used  the  source  scattering-matrix  framework  to  investigate  dipolar 
traveling  waves  that  can  be  supported  on  infinite  linear  periodic  arrays  of  lossless  penetrable 
(magnetodielectric)  spheres.  The  report  focuses  on  obtaining  the  kd  —  Pd  diagrams  for  these 
traveling  waves. 

In  future  work  we  will  attempt  to  extend  the  analysis  of  linear  periodic  arra3rs  of  pene¬ 
trable  spheres  given  in  this  report  to  two-  and  three-dimensional  periodic  arrays  of  lossless 
penetrable  spheres.  Such  an  extention  is  highly  worthwhile  in  view  of  the  facts  that  a  doubly 
negative  (DNG)  medium  can  be  formed  by  embedding  a  3D  array  of  spherical  particles  in 
a  background  medium  and  that  backward  waves  have  been  shown  in  this  report  to  be  sup¬ 
ported  by  linear  arrays  of  spheres  with  appropriately  chosen  permittivity  and  permeability. 
The  st8ui;ing  point  for  this  extension  is  equations  (32a,  32b),  and  equations  (33a,  33b)  which 
give  the  scattered  electric  and  magnetic  dipole  fields  resulting  from  an  x  directed  electric 
field  and  a  y  magnetic  field  incident  on  a  penetrable  sphere.  These  equations  then  lead  to 
equations  analogous  to  (37)  and  (39).  Instead  of  the  linear  sums  in  (37)  and  (39),  how¬ 
ever,  the  analogous  2D  and  3D  array  equations  will  contain  double  and  triple  summations 
representing  the  contributions  of  the  scattered  electric  and  magnetic  fields  of  the  spheres 
in  the  2D  and  3D  arrays  to  the  electric  and  magnetic  fields  incident  on  a  given  sphere  in 
the  array.  It  remains  to  be  seen  whether  the  complications  resulting  from  these  double  and 
triple  summations  can  be  successfully  treated. 


Figure  16:  Next-tolowest  branch  of  the  kd  —  /3d  diagram  for  a  traveling 
linear  periodic  array  of  spheres  with  Cr  =  20,  Hr  =  20,  and  ajd  =  0.45. 


Figure  17:  kd  —  I3d  diagram  for  a  traveling  wave  on  an  infinite  linear  periodic  array  of 
spheres  with  =  10, /ir  =  1>  aJid  a/d  —  0.4. 
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Figure  18:  kd  —  /3d  diagram  for  a  traveling  wave 
spheres  with  =  10,  Hr  =  1>  a/d  =  0.3. 


on  an  infinite  linear  periodic  array  of 


A  VECTOR  SPHERICAL  WAVE  FUNCTIONS 


In  this  Appendix  we  give  the  definitions  of  the  vector  spherical  wave  functions  M  and  N  used 
in  this  report.  The  definitions  are  those  used  by  Billy  Brock  [7]  adapted  for  the  exp(— iwt) 
time  dependence  used  here: 


M£(r)  =  i 


A 


(2/  +  !)(/  —  m)! 


47r/(/  +  l)(/  +  m)! 


-^zl^(kr)Pr{cos  e)eH 

sm9 


+  sin 


dX  lx=cos^ 


(64) 


NS(r)  =  i, 


(21  +  l)(/-m)! 

^  4Trl{l  +  1)(1  +  m)!  1 


+  l)J^’"(cos^)e‘"^f 
kr 


1  d 


lx=cos^ 


(65) 


I  K’H 

In  (64)  and  (65)  *md  zf^^  are  the  spherical  Bessel  and  Hankel  functions  ji  and  hi^\ 
respectively,  and  P/”(a:)  is  the  associated  Legendre  function  given  by 


Jtn 

Pr(s)  =  (-1)”(1  -a^)’”/“^P,(x),m  >  0, 

Pf'-kx)  =  (-l)”'|^^ir(:^)."‘  >  0, 


(66a) 

(66b) 


where  Pi{x)  is  the  Legendre  function 


The  vector  spherical  wave  functions  M  and  N  given  by  (64)  and  (65)  can  be  defined  in  terms 
of  the  normalized  radially-independent  vector  spherical  harmonic  function  Xjm  of  Jackson 
[6]  by 

(68) 

(69) 


M\i{T)  =  zl^\kr)XU0,(l>) 


Ji) 


and 


N£(r)  =  jVx4‘>(fa')X, „(».«, 


and  axe  related  to  the  and  vector  spherical  wave  functions  of  Hansen  [15]  by 

M,«(r)  =  iFa(r)  (70) 

and 

(71) 


NS(r)  =  iFg:,(r) 


where  the  superscript  i  =  1  when  c  =  1  and  the  superscript  i  =  2  when  c  =  3.^ 


^These  relations  between  the  vector  spherical  harmonics  of  Brock  and  those  of  Hansen  are  given  incorrectly 
in  the  original  form  of  [7]  but  have  been  corrected  by  Brock  in  an  errata. 
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B  SUMMATIONS  OF  TRIGONOMETRIC  SERIES 

In  this  Appendix  we  discuss  the  approximations  used  for  the  sums  of  the  trigonometric  series 


(72a) 

j=i  J 

(72b) 

Closed  form  expressions  are  not  available  for  these  sums.  The  IMSL  least-squares  approxi¬ 
mation  program  FNLSQ  was  used  to  compute  the  approximations 

F(a)  »  -0.1381  sin o  +  0.03212  sin 2a  -  0.9653a ln(a/7r),  0  <  a  <  tt  (73a) 


and 

G(a)  «  1.3328  -  0.1424  cos  a  +  0.01094  cos  2a  +  0.4902a^  ln(a/7r)  -  0.2417a^,  0  <  a  <  tt. 

'  (73b) 

Figures  19  and  20  show  F(a)  and  G(a),  respectively,  calculated  with  1000  terms,  along  with 
their  least  squares  approximations.  The  agreement  of  the  approximate  with  the  exact  curves 
is  excellent.  It  will  be  noted  that  although  G'(a)  =  -F(a),  when  the  approximation  for 
G(a)  is  differentiated  the  result  differs  slightly  from  the  negative  of  the  derivative  of  the 
approximation  for  F(a).  While  it  is  possible  to  obtain  an  approximation  for  F(a)  by  taking 
the  negative  of  the  derivative  of  the  approximation  for  G(a),  the  direct  least-squares  fit 
approximation  for  F(a)  that  we  have  used  distributes  the  errors  in  the  approximation  more 
uniformly  over  the  interval  from  0  to  tt. 
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